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1. Introduction 

The scenario of electrowealc baryogenesis (EWB) |jl|] owes its appeal to ttie fact that — 
unlilce most alternatives — it involves particles and interactions known to exist, and may be 
experimentally testable. (For a relatively recent review, see [||.) Baryogenesis at the weak 
scale is made possible by the fact that in the symmetric phase at temperature T, baryon 
violation is expected to occur at a rate F ~ ct^T, where is the weak fine structure 
constant [Q. Such theories require a first order weak phase transition in order to provide 
the departure from thermal equilibrium necessary for baryogenesis [Q. During the phase 
transition, bubbles of broken SUl{2) x UyiX) phase nucleate and rapidly expand, and 
the source of nonequilibrium physics necessary for baryogenesis is the interaction between 
particles in the plasma and the expanding bubble wall that separates the phases. The 
greatest challenge for an EWB mechanism is reconciling the fact that nonequilibrium CP- 
violating physics is occurring within the bubble wall itself, while baryon violation is only 
rapid outside the bubble in the symmetric phase (and perhaps in the leading edge of the 
wall). This difficulty is easily surmounted by the charge transport mechanism [|5|-P] in the 
limit when the bubble wall is narrow compared to particle mean free paths. With this 
mechanism CP violation in the phase boundary leads to the creation of net global charges 
among the particles reffecting off the domain wall; the charge is then transported into the 
region in front of the advancing domain wall. The global charge distributions in front of the 
wall contain a nonzero density of SUl{2) doublet fermions, and therefore bias anomalous 
baryon number ffuctuations toward production of baryon number 0. 

In the opposite limit, when the width of the bubble wall is large compared to the mean 
free path of fermions in the plasma, the fermions see an adiabatically changing background 
Higgs field and are not refiected. Thus the charge transport regime does not apply, and 
alternative analyses have been proposed to explain how EWB would proceed within the 
wall itself. These analyses are based on the idea of "spontaneous baryogenesis", where 
a time varying background CP-odd field arising from spontaneous symmetry breakdown 
provides the necessary CPT violation to bias B violating interactions toward the creation 

^ An alternative mechanism is proposed in [^], in which nonzero Higgs winding number is 
produced at the bubble wall, subsequently decaying in its wake, producing baryon number. 
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of net baryon number |1C]. For EWB it is a nonzero CP violating phase in a Higgs field 
that plays the role of the background field and leads to a local excess of SUl{2) doublet 
anti-fermions, causing anomalous electroweak processes to be biased in favor of producing 



baryon number within the bubble wall ||Tl|-|T3[1 . The problem with this scenario is that as 



the Higgs field turns on, baryon violating interactions turn off exponentially fast, the rate 
being proportional to exp {—2M\y/T). Taking the value of the Higgs vev where the baryon 
violation "cuts off" to be (pco, it has been argued that one can take (pco to be as large as 
(pco — ^^oi^T/g, where g is the weak coupling constant W^, in which case the resultant 



baryon asymmetry is proportional to a factor of a% ^10 ^ intrinsic to anomalous baryon 



violation. However, Dine and Thomas have argued |]T5| that in fact (j)co — oiwT/g^ and 
that as a consequence the resultant baryon asymmetry is suppressed by an additional 
factor of in the two-Higgs model, or in the minimal supersymmetric model [^| — 



suppression factors too small to account for the observed asymmetry n^/s ~ 10~^*^ when 
one includes small CP violating angles. The correct value for (^co is currently not known, 
and so the Dine-Thomas objection throws in doubt the scenario of EWB in the spontaneous 
baryogenesis (thick wall) regime. This is disturbing, since models of the electroweak sector 
that do not involve SUl{2) singlet fields generically produce wide domain walls during a 
first order phase transition, and therefore rely on the spontaneous baryogenesis mechanism 
1,0 • 



In this paper we consider the role of particle diffusion for EWB in the spontaneous 
baryogenesis regime. The potential importance of diffusion was first pointed out by Joyce, 
Prokopec and Turok [jl^ , who concluded that diffusion was inimical to spontaneous baryo- 
genesis. In the present work we arrive at a different conclusion — that diffusion moves the 
asymmetry in SUl{2) doublet fermions from within the bubble wall into the symmetric 
phase in front, where baryon violation is not exponentially suppressed. As a result the 
Dine-Thomas suppression factors are evaded, and the spontaneous baryogenesis scenario 
with diffusion ends up looking remarkably like the charge transport scenario [^p. 

We show that not only can spontaneous electroweak baryogenesis account for the 
observed baryon asymmetry of the universe, but that the asymmetry produced is not very 
sensitive to the value one takes for (pco, or to features of the weak phase transition, such 



as the width of the bubble wall (assuming it to be sufficiently wide), or its velocity. 

By analyzing the dynamical equations associated with spontaneous baryogenesis we 
are also able to treat a second important source of suppression recently discussed by 
Giudice and Shaposhnikov [|18[. They pointed out that when all of the quarks are treated 
as massless, QCD sphaleron processes equilibrate to zero the same charge densities that 
serve to drive electroweak baryogenesis. We find evidence that QCD and electroweak 
sphalerons compete, but that even in the massless quark approximation a nonzero baryon 
density results, which is surprisingly insensitive to the overall rate of sphaleron processes, 
being proportional to the ratio of electroweak to strong sphaleron rates. If this effect 
persists in a more sophisticated treatment including finite quark mass effects then one of 
the greatest uncertainties about electroweak baryogenesis — the overall rate of sphaleron 
processes — is eliminated. 

2. Diffusion equations in the two Higgs doublet model 

To be specific, we will focus on the two Higgs doublet model with non-Kobayashi- 
Maskawa CP violation as a model for electroweak baryogenesis 0, for which the sponta- 
neous baryogenesis mechanism is described in [0. During a first order phase transition 



bubbles of the broken phase will nucleate and expand. As a top quark in the plasma 
traverses the bubble wall, it interacts with the background Higgs field through its Yukawa 
coupling Ai ~ 1. In the adiabatic limit of a slowly moving or broad domain wall, the top 
quark has a spacetime dependent mass term of the form mt = XtH(r,t)e~'^^^^'^\ where 
the spacetime dependent phase 9 arises from CP violation in the Higgs potential. The 
effects of 9 are most easily analyzed by performing a hypercharge rotation of all the fields 
in the theory, removing the 9 dependence from the top quark mass term, at the expense 
of inducing an interaction of the form 

-2d^9{r,t)J^{r,t) . (2.1) 

With nonzero = —29, we see that this interaction resembles a thermodynamic charge 
constraint, with fj, being the chemical potential. Since hypercharge is not conserved in- 
side the domain wall, we refer to the dynamically generated /i as a charge potential, to 



distinguish it from an imposed constraint on a conserved charge. Exploiting the resem- 
blance, we know that the free energy is minimized by having nonzero particle densities 
ni = kiUi^T^ /6, where yi is the hypercharge, and ki a statistical factor that differs for 
fermions and bosons. In the limit that all fermions are massless except the top quark, the 
interaction (|2.1| ) leads to production of q = {tL.hL), and H particles; if fx is negative, 
then the excess of q doublets (with Y = +1/6) is also negative, and anomalous electroweak 
processes are biased to produce baryon number i. 

To see whether diffusion can be a significant effect, it is useful to consider the dimen- 
sionless quantity 

eo = D/L^Vui , 

where D is the diffusion constant for quarks, is the width of the domain wall, and 
is the wall velocity. None of these quantities is known accurately, with estimates for 
ranging from 0.1c to c/-\/3 ~ 0.6c; ranging from 10/T to 40/T, and D (classically equal 
to £c/3, where i is the mean free path) ranging from 1/T to 5/T. Evidently eo can be 
(9(1), and hence one may expect significant {i.e. 0{1)) redistributions of particle densities 
due to diffusion. As we will show, this eliminates the sensitive dependence on the value of 



discussed by Dine and Thomas [15 



In order to understand the detailed effect of the interaction ( p.l| ), we first derive the 
diffusion equation relevant for a single particle species with local number density = n 
and number current J interacting with a slowly varying background field Aq and A through 
the interaction 

C^nt = A^j^' . (2.2) 

We assume that is not exactly conserved, since is to eventually play the role of the 
hypercharge current, which is violated in the broken phase. Following the spirit of the usual 
derivation of Pick's law and the diffusion equation, we perform a simultaneous expansion 
to first order in the deviations from equilibrium (n — no), in derivatives of (n — no) and of 
the background field A, and in the particle number violating interactions. Keeping only 



Hypercharge screening proves not to affect the result [|T9|j7||; we will return to this point later. 



4 



leading terms, we find 

n = -V-J[An]-± — 

1 OLi ^2.3) 

= -v-j[An]-r ("-?°[-^]^ , 

where F is the free energy, 7 is the violation rate per unit volume of the charge Q = J d^x n, 
and we have defined the rate F as 

r - ^ • (2.4) 



The parameter k appearing in (|2.3|) is a statistical factor, 

k = (number of spin degrees of freedom) x | ^ bosons (2.5) 

I 1 fermions 

up to small corrections due to thermal particle masses. It remains to find the constituitive 
relations for J and uq. Using rotational covariance we find 

J = -DVn + aiA + a2dtA + a^VAo + ... 

(2.6) 

no = biAo + ... , 

where the ellipses refer to higher powers of derivatives acting on n and A, as well as terms 
involving the particle number violating interactions. The coefficients a and b are assumed 
to roughly equal the appropriate powers of the mean free path i dictated by dimensional 
analysis, so that the derivative expansion is justified for hydrodynamic modes of wavelength 
much longer than i. (The size of the modes we wish to consider are set by the bubble wall 
width Lw assumed to satisfy ^ £, which is the spontaneous baryogenesis regime.) 

If the current is conserved and A^ is a total divergence, then the interaction ( |2.2D has 
no physical effect. From this we deduce that 

ai = , as = —02 . 

It does not follow that hi = 0, since this coefficient only appears in eq. (|2.3| ) proportional 
to the current nonconserving effect F. We can fix 61, however, by considering the simple 
case Aq = fj, and A = 0. The minimum of the free energy is then at n = kfiT'^/G, and so 



that must be the value no in eq. (|2.3|) . It follows that the diffusion equation relevant for 
the interaction (|2.2| ) is 

n = I)V^(n-<^)-F(^-^°) 



k 



where 



no = /cAoTVe . 



This result can be readily taken over to the interaction of interest in eq. ( |2.1| ), and 
generalized to many particle species. Since = —2d^9 it follows that cp in the above 
equation vanishes, and that 



where Ui is the density of particle species i, yi is its hypercharge, and kj its statistical 
factor. The diffusion matrix Dij has diagonal elements whose magnitudes are given by the 
particle mean free paths, while the off-diagonal elements are smaller by at least one power 
of a, the fine structure constant of the transition interaction. We will work to leading 
order in perturbation theory, only keeping the diagonal entries: 



As for the r^j — hypercharge is only broken by the vev of the Higgs field, and so in the limit 
that we ignore weak mixing angles, all Yukawa couplings except for the top quark's, and 
anomalous strong interactions, the only species of particles whose densities are affected by 



handed top quark t = t^, and the Higgs particles h = {h^ +h1 + h2 + h?,). The individual 
particle numbers of these species can change through the top quark Yukawa interaction, 
the top quark mass, the Higgs self interactions, and anomalous weak interactions, at the 
rates F^, F^, Th and T^g respectively. 

Including strong sphalerons (with a rate Tgs) would allow the generation of right 
handed bottom quarks, b = bn, as well as first and second family quarks. However since 
strong sphalerons are the only processes which generate significant numbers of first and 
second family quarks, and all quarks have approximately the same diffusion constant, these 
densities are constrained algebraically in terms of b to satisfy 




(2.7) 



= DAj + 0{a) . 



a nonzero 9 are the left handed third family doublet denoted by = (t^ + 6l), the right 



qiL = <?2L = -2wi? = -2dR = -2sr = -2cr = -2bR = -2b 
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(2.8) 



To simplify the equations we may use the fact that the only charge potential generated 
is for hypercharge, and thus each no is proportional to k times the hypercharge for the 
species — this allows us to eliminate all no in favor of ho. Using the relations 



^qo/kq + lo/ki oc [Sy^ - yi) = (1/6)3 -1/2 = 

{qo/kq - ho/kh - to/kt) (X {y^ - yn - yt) = (1/6) - (-1/2) - (2/3) = 

we see that the source terms proportional to 6 drop out of Vy, Vyjs interactions and are 
only present proportional to the hypercharge violating and T^, consistent with the 
arguments of Dine and Thomas . 
Eq. ( p.7|) may now be written as 



q = DqV'^q - Ty [q/kq - h/ku - t/h] - [q/kq - t/kt - /lo/^h 

- 6r^, [3(q - Ah)/kq + l/h] - 6r,, [2q/kq - t/h - %/ki,] 
i = DtVH - Ty [-q/kq + h/kh + t/kt] 

- Vm [-q/kq + t/kt + ho/kh] + 3r,, [2q/kq - t/h - 9b/kt 
h = DhV^h - Ty [-q/kq + t/kt + h/kh] - Th{h - ho)/kh 

b = DtV% + 3Tss [2q/kq - t/h - %/kb] 

where the k factors are given by (|2.5| ) 



(2.9) 



kq = 6, A;i = 3, kh = ^. kb = 3 . (2.10) 



Several simplifications of equations ( |2.9|) can be made. First we ignore the curvature 
of the bubble wall, and so F^, F/i, and F^g are only functions of z = \r — Vwt\, where 
Vw is the bubble wall velocity. Explicitly, we assume that F^, and F^ vary as the square 
of the Higgs vev, while F^^s is approximated by a step function: r^^s{z) — 9{z — Zco)Tws- 
The parameters Dt, Dh, Fy and Tss are taken to be spacetime constants. We will assume 
that the density perturbations of interest are functions of z and t only. We will also take 
Dq Dt — -Db, which is correct up to order aw/o^s. Then, ignoring (9(F^g) effects on the 
final baryon density, we may set b = —{q + t) and eliminate the b equation; furthermore, 
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since no leptons are produced in the limit that we ignore lepton Yukawa couphngs we may 
set / = 0. 

With these assumptions we arrive at the diffusion equations for q{z,t),t{z,t), and 
h{z, t) in the rest frame of the bubble wall (where z is the coordinate normal to the 
surface): 

q = v^d.q + DqV^q - Ty [q/6 - h/8 - t/3] - [q/6 - t/3 - ho/8] 

- 3T^s [5q + 4t] - 2Tss [Wq + 8t] 
i = v^d.t + Dqdlq + Vy [q/Q - h/8 - t/3] (2.11) 

+ [q/Q - t/3 - ho/8] + r,, [lOq + 8t] 

h = v^d.h + Dhdlh + +Ty [q/% - h/8 - t/3] - Th{z) [{h - ho)/8] . 
In our past treatment of this system [|l^ we assumed that h and q — t (axial top 
number) had time independent solutions in the wall rest frame equal to the values which 
minimize the free energy with zero baryon number. This was tantamount to treating the 
dimensionless numbers Vzco/vw and Tz'^^/D as being much larger than unity. If Zco is very 
small, as suggested by Dine and Thomas [l^], than this assumption is false — especially 
since the rates and Fh are proportional to the Higgs vev and therefore vanish at the 
leading edge of the wall. Eq. ( p.ll|) allows us to deal with realistic values for F, D, and 



3. Solution of the rate equations 

In this paper we will only use reasonable estimates of the diffusion constants D and 
rates T{z) in order to demonstrate the general properties of the solutions to Eq. ( |2.11| ). 
We will take the bubble wall profile at a transition temperature T to have the form 

{H{z)) ~ Te-'^^^\l - tanh{z/L^))/2 , 

(3.1) 

e{z) = Ae{l - tanh(2/L^))/2 , 



where AO is a measure of the size of CP violation in the theory. It follows from (|2.7| ) that 
the charge potential ho in eq. ( |2.11|) is given by 

ho = -Ti'^wT — (3.2) 
3 az 



and that the solutions wiU scale linearly with A^. For the diffusion constants we take 



Dt = D, = 3/T , Dh = 10/T , (3.3) 

where we have used the estimates from for Dq and have scaled Dh — Dq{as/ cxw)- For 
the rates (defined in Eq. (|2.4|) as F = 67/T"^) we take 

F;, ~ F^ = {Hiz))^X^jT , (3.4) 

and 

F^s = 6«;q;^T , F^^ = Gn-a'^T . (3.5) 

The steady state solution to Eqs. ( |2.11| ) was found by numerical integration. The 
results for the choice of parameters k = 1, = 0.1, = 10/T and A9 = —tv are given 
in fig. ^ The results clearly indicate that diffusion plays a significant role, since nonzero 
particle densities are seen to extend well beyond the bubble wall into the symmetric phase. 
An effect of strong sphalerons is seen in the nonzero q+t density, which vanishes as F^g — > 0. 
Solutions to ( p.ll|) were also obtained for a variety of wall velocities (0.05 < < 0.3), 
wall widths {10/T < < 40/T), and sphaleron rates (0.1 < k < 100). 

The baryon density deep within the bubble {z large and negative in these coordinates) 
is computed by simply integrating the equation for total lefthanded fermion number (equal 
to (5^ + 4t), after solving algebraically for the number densities of light quarks) from far 
outside the bubble down to Zco- 

3F 

UB^^^l dz {f>q{z) + U{z)) . (3.6) 



As discussed in the introduction, the value of Zco is controversial, and a large value (cor- 
responding to small (j)co) suppresses baryon number when diffusion is ignored. In fig. || 
we present the final baryon asymetry ub /s of the universe as a function of Zco for several 
values of k, as well as the solution for k = 1 when all the diffusion constants are taken to 
be very small, Dq = = 0.1/T. Apparently spontaneous electroweak baryogenesis can 
quite easily account for the observed asymmetry ns/s ~ lO"^'^ with CP violation at the 
~ 10-2 - 10-3 level. 
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There are several remarkable features to fig. ^: First, it illustrates the Dine-Thomas 
sensitivity to the value of Zco when diffusion is ignored, but the sensitivity is practically 
eliminated when diffusion is taken into account. A second amusing feature is that the final 
baryon asymmetry turns out to be approximately independent of k in the range examined. 
Ignoring the effects of strong sphalerons one would have found oc k, but as pointed 
out by Giudice and Shaposhnikov [^| the left handed baryon number is driven to zero 



by strong sphalerons in the limit that the thermal fermion masses are neglected. We find 
numerically that the baryon number with Zco = 3.0L^(, is inversely proportional to Tgs for 
0.5 < K < 100, cancelling the linear k dependence for electroweak sphalerons. For k < 0.1, 
we find ub to be insensitive to Tss and proportional to k (an intuitive understanding of this 
number is that strong sphalerons are irrelevant when the timescale l/Tss is much longer 
than the typical time D/v'^ that the diffusing fermions spend in the symmetric phase.) i. 



We conclude this section by noting that solutions to Eq. (|2.11|) with k = 1 and a 
variety of wall velocities and wall widths reveal that the baryon asymmetry one creates 
is not extremely sensitive to either or Lw We find that for Lyj = 10/T and 
ranging from 0.05 to 0.3, the baryon asymmetry changes by no more than a factor of 
2; for a broader wall with v^^ = 0.1 and Lyj = 40/T, ub was decreased by a factor of 
~ 5. When we artificially set the diffusion constants and Tss to be very small and assume 
(j)co = i4:CtwT/g, K = 1, we reproduce our earlier results which neglected both diffusion 
and strong sphalerons 0. 

4. Outlook 

We conclude that a dynamical treatment of spontaneous electroweak baryogenesis with 
finite interaction rates, nonzero diffusion constants, and strong sphaleron effects eliminates 



^ We have not presented data for k > 5.0 since the neglected mt and a corrections to the 
relation n ~ kfiT'^/d discussed in are expected to become more important as k gets larger. 



These effects are expected to further enhance the baryon asymmetry. Note however that even 
neglecting these effects one finds a sizable baryon asymmetry in our dynamical calculation, unlike 



in the equilibrium calculation of |18] 



^ Our review contains a typographical error in Eq. 37: the right side of Eq. 37 should 
include the factor (A/"/. 33) instead of (A/"/.!) 
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or moderates large suppression factors found in an equilibrium calculation. Our work 
indicates that the EWB scenario remains a viable explanation for the baryon asymmetry 
of the universe. Furthermore, the results appear not to be very sensitive to the bubble 
wall velocity or width, to the point in the wall where anomalous baryon violation becomes 
suppressed, or to the sphaleron rate parameter k for a broad range of these parameters. 
This is encouraging, since these quantities are poorly known even in specific models. On 
the other hand, our results are rather sensitive to the values of the diffusion lengths and 
interaction rates, which we have only estimated. These quantities are computable from 
known physics, and this should be done before making firm conclusions. Corrections of 
order and a to the relationship between density and charge potential should also be 
included, particularly for large k. 

A more complete treatment would also include the effects of hypercharge screening 
1^,^. Screening is known not to affect the spontaneous baryogenesis in equilibrium cal- 



culations, since three quark doublets are drawn in for every anti-lepton doublet, so that 
there is no net effect on weak sphalerons. In a dynamical calculation one might expect the 
anti-leptons to screen more efficiently since they have a longer mean free path, in which 
case the baryon number produced would only be enhanced. 

Finally, we note that diffusion will be even more important for spontaneous baryoge- 



nesis in the minimal supersymmetric model JI^, due to the contributions from charginos 
and neutralinos which have much longer mean free paths than quarks. In addition, strong 
sphaleron suppression will be absent due the combinations of charges that bias baryon 
number in that model. 
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Figure Captions 



Fig. 1. Time independent particle number densities normalized to the entropy as a func- 
tion of position in the rest frame of the expanding bubble wall. The parameters 
used were Vw = 0.1, = 10/T, k = 1, and A9 = — tt; the D and F coefficients 
used are given in Eqs. (|3.3|) , (|3^), ( p.5|) . The densities /i, q — t, and q + t refer to 
Higgs number, tL + Bl — ti?, and + 6l + ^i? respectively. The midpoint of the 
bubble wall is at ^ = 0. 

Fig. 2. The final baryon to entropy ratio of the universe for several values of k, plotted 
as LogiQ[{n B/s){g^/ 100) {— TV / A9)] versus Zco, the point where weak sphalerons 
become exponentially suppressed. Zco = —O-bL^j corresponds to (j)co — 14:awT/g, 
while Zco = 3.0L^y corresponds to (^co — o-wT/g. The dashed line is the solution to 
Eq. ( ]2.11D assuming very small diffusion constants, illustrating the Dine-Thomas 



effect. 
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